Novel regulators of islet function identified from genetic variation in mouse islet Ca2+ oscillations.

Insufficient insulin secretion to meet metabolic demand results in diabetes. The intracellular flux of Ca2+ into ß-cells triggers insulin release. Since genetics strongly influences variation in islet secretory responses, we surveyed islet Ca2+ dynamics in eight genetically diverse mouse strains. We found high strain variation in response to four conditions: (1) 8 mM glucose; (2) 8 mM glucose plus amino acids; (3) 8 mM glucose, amino acids, plus 10 nM glucose-dependent insulinotropic polypeptide (GIP); and (4) 2 mM glucose. These stimuli interrogate ß-cell function, α- to ß-cell signaling, and incretin responses. We then correlated components of the Ca2+ waveforms to islet protein abundances in the same strains used for the Ca2+ measurements. To focus on proteins relevant to human islet function, we identified human orthologues of correlated mouse proteins that are proximal to glycemic-associated single-nucleotide polymorphisms in human genome-wide association studies. Several orthologues have previously been shown to regulate insulin secretion (e.g. ABCC8, PCSK1, and GCK), supporting our mouse-to-human integration as a discovery platform. By integrating these data, we nominate novel regulators of islet Ca2+ oscillations and insulin secretion with potential relevance for human islet function. We also provide a resource for identifying appropriate mouse strains in which to study these regulators.

An Enhancer Within Abcb11 Regulates G6pc2 in C57BL/6 Mouse Pancreatic Islets.

G6PC2 is predominantly expressed in pancreatic islet ß-cells where it encodes a glucose-6-phosphatase catalytic subunit that modulates the sensitivity of insulin secretion to glucose by opposing the action of glucokinase, thereby regulating fasting blood glucose (FBG). Prior studies have shown that the G6pc2 promoter alone is unable to confer sustained islet-specific gene expression in mice, suggesting the existence of distal enhancers that regulate G6pc2 expression. Using information from both mice and humans and knowledge that single nucleotide polymorphisms (SNPs) both within and near G6PC2 are associated with variations in FBG in humans, we identified several putative enhancers 3' of G6pc2. One region, herein referred to as enhancer I, resides in the 25th intron of Abcb11 and binds multiple islet-enriched transcription factors. CRISPR-mediated deletion of enhancer I in C57BL/6 mice had selective effects on the expression of genes near the G6pc2 locus. In isolated islets, G6pc2 and Spc25 expression were reduced ∼50%, and Gm13613 expression was abolished, whereas Cers6 and nostrin expression were unaffected. This partial reduction in G6pc2 expression enhanced islet insulin secretion at basal glucose concentrations but did not affect FBG or glucose tolerance in vivo, consistent with the absence of a phenotype in G6pc2 heterozygous C57BL/6 mice.

Lipidomic QTL in Diversity Outbred mice identifies a novel function for α/ß hydrolase domain 2 (Abhd2) as an enzyme that metabolizes phosphatidylcholine and cardiolipin.

We and others have previously shown that genetic association can be used to make causal connections between gene loci and small molecules measured by mass spectrometry in the bloodstream and in tissues. We identified a locus on mouse chromosome 7 where several phospholipids in liver showed strong genetic association to distinct gene loci. In this study, we integrated gene expression data with genetic association data to identify a single gene at the chromosome 7 locus as the driver of the phospholipid phenotypes. The gene encodes α/ß-hydrolase domain 2 (Abhd2), one of 23 members of the ABHD gene family. We validated this observation by measuring lipids in a mouse with a whole-body deletion of Abhd2. The Abhd2KO mice had a significant increase in liver levels of phosphatidylcholine and phosphatidylethanolamine. Unexpectedly, we also found a decrease in two key mitochondrial lipids, cardiolipin and phosphatidylglycerol, in male Abhd2KO mice. These data suggest that Abhd2 plays a role in the synthesis, turnover, or remodeling of liver phospholipids.

Lipidomic QTL in Diversity Outbred mice identifies a novel function for α/ß hydrolase domain 2 ( Abhd2 ) as an enzyme that metabolizes phosphatidylcholine and cardiolipin.

We and others have previously shown that genetic association can be used to make causal connections between gene loci and small molecules measured by mass spectrometry in the bloodstream and in tissues. We identified a locus on mouse chromosome 7 where several phospholipids in liver showed strong genetic association to distinct gene loci. In this study, we integrated gene expression data with genetic association data to identify a single gene at the chromosome 7 locus as the driver of the phospholipid phenotypes. The gene encodes α/ß-hydrolase domain 2 ( Abhd2 ), one of 23 members of the ABHD gene family. We validated this observation by measuring lipids in a mouse with a whole-body deletion of Abhd2 . The Abhd2 KO mice had a significant increase in liver levels of phosphatidylcholine and phosphatidylethanolamine. Unexpectedly, we also found a decrease in two key mitochondrial lipids, cardiolipin and phosphatidylglycerol, in male Abhd2 KO mice. These data suggest that Abhd2 plays a role in the synthesis, turnover, or remodeling of liver phospholipids.

ß Cell-specific deletion of Zfp148 improves nutrient-stimulated ß cell Ca2+ responses.

Insulin secretion from pancreatic ß cells is essential for glucose homeostasis. An insufficient response to the demand for insulin results in diabetes. We previously showed that ß cell-specific deletion of Zfp148 (ß-Zfp148KO) improves glucose tolerance and insulin secretion in mice. Here, we performed Ca2+ imaging of islets from ß­Zfp148KO and control mice fed both a chow and a Western-style diet. ß-Zfp148KO islets demonstrated improved sensitivity and sustained Ca2+ oscillations in response to elevated glucose levels. ß-Zfp148KO islets also exhibited elevated sensitivity to amino acid-induced Ca2+ influx under low glucose conditions, suggesting enhanced mitochondrial phosphoenolpyruvate-dependent (PEP-dependent), ATP-sensitive K+ channel closure, independent of glycolysis. RNA-Seq and proteomics of ß-Zfp148KO islets revealed altered levels of enzymes involved in amino acid metabolism (specifically, SLC3A2, SLC7A8, GLS, GLS2, PSPH, PHGDH, and PSAT1) and intermediary metabolism (namely, GOT1 and PCK2), consistent with altered PEP cycling. In agreement with this, ß-Zfp148KO islets displayed enhanced insulin secretion in response to l-glutamine and activation of glutamate dehydrogenase. Understanding pathways controlled by ZFP148 may provide promising strategies for improving ß cell function that are robust to the metabolic challenge imposed by a Western diet.

Coding variants identified in patients with diabetes alter PICK1 BAR domain function in insulin granule biogenesis.

Bin/amphiphysin/Rvs (BAR) domains are positively charged crescent-shaped modules that mediate curvature of negatively charged lipid membranes during remodeling processes. The BAR domain proteins PICK1, ICA69, and the arfaptins have recently been demonstrated to coordinate the budding and formation of immature secretory granules (ISGs) at the trans-Golgi network. Here, we identify 4 coding variants in the PICK1 gene from a whole-exome screening of Danish patients with diabetes that each involve a change in positively charged residues in the PICK1 BAR domain. All 4 coding variants failed to rescue insulin content in INS-1E cells upon knock down of endogenous PICK1. Moreover, 2 variants showed dominant-negative properties. In vitro assays addressing BAR domain function suggested that the coding variants compromised BAR domain function but increased the capacity to cause fission of liposomes. Live confocal microscopy and super-resolution microscopy further revealed that PICK1 resides transiently on ISGs before egress via vesicular budding events. Interestingly, this egress of PICK1 was accelerated in the coding variants. We propose that PICK1 assists in or complements the removal of excess membrane and generic membrane trafficking proteins, and possibly also insulin, from ISGs during the maturation process; and that the coding variants may cause premature budding, possibly explaining their dominant-negative function.

Identification of direct transcriptional targets of NFATC2 that promote ß cell proliferation.

The transcription factor NFATC2 induces ß cell proliferation in mouse and human islets. However, the genomic targets that mediate these effects have not been identified. We expressed active forms of Nfatc2 and Nfatc1 in human islets. By integrating changes in gene expression with genomic binding sites for NFATC2, we identified approximately 2200 transcriptional targets of NFATC2. Genes induced by NFATC2 were enriched for transcripts that regulate the cell cycle and for DNA motifs associated with the transcription factor FOXP. Islets from an endocrine-specific Foxp1, Foxp2, and Foxp4 triple-knockout mouse were less responsive to NFATC2-induced ß cell proliferation, suggesting the FOXP family works to regulate ß cell proliferation in concert with NFATC2. NFATC2 induced ß cell proliferation in both mouse and human islets, whereas NFATC1 did so only in human islets. Exploiting this species difference, we identified approximately 250 direct transcriptional targets of NFAT in human islets. This gene set enriches for cell cycle-associated transcripts and includes Nr4a1. Deletion of Nr4a1 reduced the capacity of NFATC2 to induce ß cell proliferation, suggesting that much of the effect of NFATC2 occurs through its induction of Nr4a1. Integration of noncoding RNA expression, chromatin accessibility, and NFATC2 binding sites enabled us to identify NFATC2-dependent enhancer loci that mediate ß cell proliferation.

INFIMA leverages multi-omics model organism data to identify effector genes of human GWAS variants.

Genome-wide association studies reveal many non-coding variants associated with complex traits. However, model organism studies largely remain as an untapped resource for unveiling the effector genes of non-coding variants. We develop INFIMA, Integrative Fine-Mapping, to pinpoint causal SNPs for diversity outbred (DO) mice eQTL by integrating founder mice multi-omics data including ATAC-seq, RNA-seq, footprinting, and in silico mutation analysis. We demonstrate INFIMA's superior performance compared to alternatives with human and mouse chromatin conformation capture datasets. We apply INFIMA to identify novel effector genes for GWAS variants associated with diabetes. The results of the application are available at http://www.statlab.wisc.edu/shiny/INFIMA/ .

From methylene bridged diindole to carbonyl linked benzimidazoleindole: Development of potent and metabolically stable PCSK9 modulators.

Proprotein convertase subtilisin/kexin type 9 (PCSK9) is a recently validated therapeutic target for lowering low-density lipoprotein cholesterol (LDL-C). Through phenotypic screening, we previously discovered a class of small-molecules with a 2,3'-diindolymethane (DIM) skeleton that can decrease the expression of PCSK9. But these compounds have low potency and low metabolically stability. After performing structure-activity relationship (SAR) optimization by nitrogen scan, deuterium substitution and fluorine scan, we identified a series of much more potent and metabolically stable PCSK9 modulators. A preliminary in vivo pharmacokinetic study was performed for representative analogues difluorodiindolyketone (DFDIK) 12 and difluorobenzoimidazolylindolylketone (DFBIIK-1) 13. The in vitro metabolic stability correlate well with the in vivo data. The most potent compound 21 has the EC50 of 0.15 nM. Our SAR studies also indicated that the NH on the indole ring of 21 can tolerate more function groups, which may facilitate the mechanism of action studies and also allow further improvement of the pharmacological properties.

Secretion of Recombinant Interleukin-22 by Engineered Lactobacillus reuteri Reduces Fatty Liver Disease in a Mouse Model of Diet-Induced Obesity.

The incidence of metabolic syndrome continues to rise globally. In mice, intravenous administration of interleukin-22 (IL-22) ameliorates various disease phenotypes associated with diet-induced metabolic syndrome. In patients, oral treatment is favored over intravenous treatment, but methodologies to deliver IL-22 via the oral route are nonexistent. The goal of this study was to assess to what extent engineered Lactobacillus reuteri secreting IL-22 could ameliorate nonalcoholic fatty liver disease. We used a mouse model of diet-induced obesity and assessed various markers of metabolic syndrome following treatment with L. reuteri and a recombinant derivative. Mice that received an 8-week treatment of wild-type probiotic gained less weight and had a smaller fat pad than the control group, but these phenotypes were not further enhanced by recombinant L. reuteri However, L. reuteri secreting IL-22 significantly reduced liver weight and triglycerides at levels that exceeded those of the probiotic wild-type treatment group. Our findings are interesting in light of the observed phenotypes associated with reduced nonalcoholic liver disease, in humans the most prevalent chronic liver disease, following treatment of a next-generation probiotic that is administered orally. Once biological and environmental containment strategies are in place, therapeutic applications of recombinant Lactobacillus reuteri are on the horizon.IMPORTANCE In humans, nonalcoholic fatty liver disease (NAFLD) is the most prevalent liver disease due to the increased prevalence of obesity. While treatment of NAFLD is often geared toward lifestyle changes, such as diet and exercise, the use of dietary supplements such as probiotics is underinvestigated. Here, we report that probiotic Lactobacillus reuteri reduces fatty liver in a mouse model of diet-induced obesity. This phenotype was further enhanced upon delivery of recombinant interleukin-22 by engineered Lactobacillus reuteri These observations pave the road to a better understanding of probiotic mechanisms driving the reduction of diet-induced steatosis and to development of next-generation probiotics for use in the clinic. Ultimately, these studies may lead to rational selection of (engineered) probiotics to ameliorate fatty liver disease.

Gene loci associated with insulin secretion in islets from non-diabetic mice.

Genetic susceptibility to type 2 diabetes is primarily due to ß-cell dysfunction. However, a genetic study to directly interrogate ß-cell function ex vivo has never been previously performed. We isolated 233,447 islets from 483 Diversity Outbred (DO) mice maintained on a Western-style diet, and measured insulin secretion in response to a variety of secretagogues. Insulin secretion from DO islets ranged >1,000-fold even though none of the mice were diabetic. The insulin secretory response to each secretagogue had a unique genetic architecture; some of the loci were specific for one condition, whereas others overlapped. Human loci that are syntenic to many of the insulin secretion QTL from mouse are associated with diabetes-related SNPs in human genome-wide association studies. We report on three genes, Ptpn18, Hunk and Zfp148, where the phenotype predictions from the genetic screen were fulfilled in our studies of transgenic mouse models. These three genes encode a non-receptor type protein tyrosine phosphatase, a serine/threonine protein kinase, and a Krϋppel-type zinc-finger transcription factor, respectively. Our results demonstrate that genetic variation in insulin secretion that can lead to type 2 diabetes is discoverable in non-diabetic individuals.

Genetic Drivers of Pancreatic Islet Function.

The majority of gene loci that have been associated with type 2 diabetes play a role in pancreatic islet function. To evaluate the role of islet gene expression in the etiology of diabetes, we sensitized a genetically diverse mouse population with a Western diet high in fat (45% kcal) and sucrose (34%) and carried out genome-wide association mapping of diabetes-related phenotypes. We quantified mRNA abundance in the islets and identified 18,820 expression QTL. We applied mediation analysis to identify candidate causal driver genes at loci that affect the abundance of numerous transcripts. These include two genes previously associated with monogenic diabetes (PDX1 and HNF4A), as well as three genes with nominal association with diabetes-related traits in humans (FAM83E, IL6ST, and SAT2). We grouped transcripts into gene modules and mapped regulatory loci for modules enriched with transcripts specific for α-cells, and another specific for δ-cells. However, no single module enriched for ß-cell-specific transcripts, suggesting heterogeneity of gene expression patterns within the ß-cell population. A module enriched in transcripts associated with branched-chain amino acid metabolism was the most strongly correlated with physiological traits that reflect insulin resistance. Although the mice in this study were not overtly diabetic, the analysis of pancreatic islet gene expression under dietary-induced stress enabled us to identify correlated variation in groups of genes that are functionally linked to diabetes-associated physiological traits. Our analysis suggests an expected degree of concordance between diabetes-associated loci in the mouse and those found in human populations, and demonstrates how the mouse can provide evidence to support nominal associations found in human genome-wide association mapping.

The Transcription Factor Nfatc2 Regulates ß-Cell Proliferation and Genes Associated with Type 2 Diabetes in Mouse and Human Islets.

Human genome-wide association studies (GWAS) have shown that genetic variation at >130 gene loci is associated with type 2 diabetes (T2D). We asked if the expression of the candidate T2D-associated genes within these loci is regulated by a common locus in pancreatic islets. Using an obese F2 mouse intercross segregating for T2D, we show that the expression of ~40% of the T2D-associated genes is linked to a broad region on mouse chromosome (Chr) 2. As all but 9 of these genes are not physically located on Chr 2, linkage to Chr 2 suggests a genomic factor(s) located on Chr 2 regulates their expression in trans. The transcription factor Nfatc2 is physically located on Chr 2 and its expression demonstrates cis linkage; i.e., its expression maps to itself. When conditioned on the expression of Nfatc2, linkage for the T2D-associated genes was greatly diminished, supporting Nfatc2 as a driver of their expression. Plasma insulin also showed linkage to the same broad region on Chr 2. Overexpression of a constitutively active (ca) form of Nfatc2 induced ß-cell proliferation in mouse and human islets, and transcriptionally regulated more than half of the T2D-associated genes. Overexpression of either ca-Nfatc2 or ca-Nfatc1 in mouse islets enhanced insulin secretion, whereas only ca-Nfatc2 was able to promote ß-cell proliferation, suggesting distinct molecular pathways mediating insulin secretion vs. ß-cell proliferation are regulated by NFAT. Our results suggest that many of the T2D-associated genes are downstream transcriptional targets of NFAT, and may act coordinately in a pathway through which NFAT regulates ß-cell proliferation in both mouse and human islets.

Obesity and genetics regulate microRNAs in islets, liver, and adipose of diabetic mice.

Type 2 diabetes results from severe insulin resistance coupled with a failure of b cells to compensate by secreting sufficient insulin. Multiple genetic loci are involved in the development of diabetes, although the effect of each gene on diabetes susceptibility is thought to be small. MicroRNAs (miRNAs) are noncoding 19-22-nucleotide RNA molecules that potentially regulate the expression of thousands of genes. To understand the relationship between miRNA regulation and obesity-induced diabetes, we quantitatively profiled approximately 220 miRNAs in pancreatic islets, adipose tissue, and liver from diabetes-resistant (B6) and diabetes-susceptible (BTBR) mice. More than half of the miRNAs profiled were expressed in all three tissues, with many miRNAs in each tissue showing significant changes in response to genetic obesity. Furthermore, several miRNAs in each tissue were differentially responsive to obesity in B6 versus BTBR mice, suggesting that they may be involved in the pathogenesis of diabetes. In liver there were approximately 40 miRNAs that were downregulated in response to obesity in B6 but not BTBR mice, indicating that genetic differences between the mouse strains play a critical role in miRNA regulation. In order to elucidate the genetic architecture of hepatic miRNA expression, we measured the expression of miRNAs in genetically obese F2 mice. Approximately 10% of the miRNAs measured showed significant linkage (miR-eQTLs), identifying loci that control miRNA abundance. Understanding the influence that obesity and genetics exert on the regulation of miRNA expression will reveal the role miRNAs play in the context of obesity-induced type 2 diabetes.

Combined expression trait correlations and expression quantitative trait locus mapping.

Coordinated regulation of gene expression levels across a series of experimental conditions provides valuable information about the functions of correlated transcripts. The consideration of gene expression correlation over a time or tissue dimension has proved valuable in predicting gene function. Here, we consider correlations over a genetic dimension. In addition to identifying coregulated genes, the genetic dimension also supplies us with information about the genomic locations of putative regulatory loci. We calculated correlations among approximately 45,000 expression traits derived from 60 individuals in an F2 sample segregating for obesity and diabetes. By combining the correlation results with linkage mapping information, we were able to identify regulatory networks, make functional predictions for uncharacterized genes, and characterize novel members of known pathways. We found evidence of coordinate regulation of 174 G protein-coupled receptor protein signaling pathway expression traits. Of the 174 traits, 50 had their major LOD peak within 10 cM of a locus on Chromosome 2, and 81 others had a secondary peak in this region. We also characterized a Riken cDNA clone that showed strong correlation with stearoyl-CoA desaturase 1 expression. Experimental validation confirmed that this clone is involved in the regulation of lipid metabolism. We conclude that trait correlation combined with linkage mapping can reveal regulatory networks that would otherwise be missed if we studied only mRNA traits with statistically significant linkages in this small cross. The combined analysis is more sensitive compared with linkage mapping alone.